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Abstract 

New physical effects emerge from an interplay between the electron parity number and persistent currents in superconducting 
nanorings. An odd electron, being added to the ring, produces a countercurrent which may substantially modify the ground 
state properties of the system. In superconducting nanorings with an embedded normal metal layer a novel U n/N -junction" 
state can occur for the odd number of electrons. Changing this number from even to odd yields spontaneous supercurrent 
in the ground state of such rings without any externally applied magnetic flux. Further peculiar features of the parity effect 
are expected in structures with resonant electron transport across the weak link. 
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1. Introduction 

It is well established both theoretically and exper- 
imentally [1,2,3,4,5,6] that thermodynamic properties 
of isolated superconducting systems are sensitive to 
the parity of the total number of electrons M despite 
this number may be macroscopically large. The funda- 
mental difference between superconductors with even 
and odd N is due to the fact that the number of elec- 
trons in the condensate of Cooper pairs should neces- 
sarily be even. Hence, for odd N at least one electron 
remains unpaired having an extra energy equal to the 
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superconducting gap A. This difference - hindered by 
the entropy factor at sufficiently high temperatures - 
becomes important in the low temperature limit. 

Can the supercurrent be affected by this parity ef- 
fect? Recently we have argued [7] that (a) the answer 
to this question is positive and (b) in superconducting 
nanorings interrupted by a weak link parity-restricted 
supercurrent can strongly deviate from one evaluated 
for the grand canonical ensemble. In this paper we will 
extend our previous analysis focusing it on supercon- 
ducting rings with embedded piece of a normal metal 
(SNS rings). We will consider transparent nanojunc- 
tions as well as structures with low transmission of NS- 
interfaces. In the latter case resonant effects become 
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important leading to substantial modifications of the 
physical picture of the parity effect [7] . 



2. Parity projection formalism 

In the subsequent analysis of parity-restricted per- 
sistent currents (PC) in superconducting nanorings we 
will employ the parity projection formalism [4,5,6]. Re- 
capitulating the key points of this approach we will 
closely follow Ref . [5] . We will then derive a general for- 
mula for the circulating supercurrent in isolated super- 
conducting rings with even/odd number of electrons. 

The grand canonical partition function Z(T,[i) = 
Tre -/3(H-MA0 ig linked t0 the can onical one Z(T, A/") as 



(1) 



Here and below H is the system Hamiltonian and (3 = 
1/T. Inverting this relation and defining the canoni- 
cal partition functions Z e and Z respectively for even 
(Af = Af e ) and odd (Af = Afo) ensembles, one gets 
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are the parity projected grand canonical partition func- 
tions. For Af ^> 1 it is sufficient to evaluate the integral 
in (2) within the saddle point approximation 
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are the parity projected thermodynamic potentials, 
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and 
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"Chemical potentials" /j, e / a are defined by the saddle 
point condition 
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-dQ e / (T, n e / )/du e / . 



Let us emphasize that fi/ is just the standard grand 
canonical thermodynamic potential and fit represents 
the corresponding potential linked to the auxiliary par- 
tition function Z(T, fi+inT). This function can be con- 
veniently evaluated by finding the true grand canon- 
ical partition function Z(T,[i), expressing the result 
as a sum over the Fermi Matsubara frequencies uif — 
2nT(l + 1/2) and then substituting the Bose Matsub- 
ara frequencies cut = 27rT/ (I — 0, ±1, ...) instead of LOf. 

Consider now an isolated superconducting ring 
pierced by the magnetic flux $. Making use of the 
above expressions one one can easily determine the 
equilibrium current 7 e / circulating inside the ring: 
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where the upper/lower sign corresponds to the 
even/odd ensemble and we have defined 
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Below we will make use of these equations in order to 
evaluate PC in superconducting nanorings. 



3. Superconducting rings with nanojunctions 

3.1. The model and general analysis 

Before we turn to concrete calculations let us specify 
the model for our system. In what follows we will in- 
vestigate PC in superconducting nanorings with cross 
section s and perimeter L — 2nR. Superconducting 
properties of such rings will be described within the 
(parity projected) mean field BCS theory. A necessary 
- though not always sufficient - validity condition for 
this mean field approach reads N T ~ p%s 3> 1, i.e. 
the total number of conducting channels in the ring N T 
should remain large. More detailed requirements can 
be formulated with the aid of the results [10,11]. For 
generic wires QPS effects are small provided the pa- 
rameter yfs exceeds ~ 10 nm [10,11,12,13]. Consider- 
ing then homogeneous rings one can prove [7] that in 
such rings the effect of the electron parity number on 
PC remains small in the parameter ~ 1/N T <C 1. On 
the other hand, in rings with only few number of chan- 
nels [8,9] iV r ~ 1 fluctuation effects may entirely wipe 
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Fig. 1. Superconducting ring with embedded SNS junction 
of length d. 

out the supercurrent thus making the mean field ap- 
proach obsolete. 

The way out is to consider superconducting rings 
with JV r 3> 1 interrupted by a weak link with only few 
number of conducting channels N. In such systems the 
mean field BCS description remains applicable and, on 
the other hand, the parity effect can be large due to 
the condition N ~ 1. Quantum point contacts (QPC) 
or, more generally, SNS junctions (of length d) can be 
used for practical realization of such weak links. The 
corresponding structure is depicted in Fig. 1. 

In this case the thermodynamic potential of the sys- 
tem Q consists of two different contributions [141 



n = n (r) (fj,, t, $ x>v >) + n {cl ( M , t, <p) 



(c), 



(7) 



respectively from the bulk part of the ring and from the 
contact. The optimal value of the phase difference ip 
across the weak link is fixed by the condition dtt / dp = 
which reads 



2e dQ (c) 
h dip 



(8) 



In (8) we made use of the fact that the thermodynamic 
potentials of the ring depend both on & x and ip only via 
the superfiuid velocity vs = (l/AimiR)((p— 27r$ I /<E>o), 
in which case one can put d/d& x — — (2e/hc)(d/dip). 
Here and below <E>o is the superconducting flux quan- 
tum. The left-hand side of Eq. (8) represents the cur- 
rent flowing inside the superconducting ring 1^ = 
-cdn (r) d$> x ~ {ev F N v /L)(p - 27r$ ;c /$o). This value 
should be equal to the current across the weak link 
which is given by the right-hand side of Eq. (8). Es- 
timating the maximum value of the latter for a single 
channel contact as 2eTA/h, we obtain 



p~2n^, if 

<p ~ 27m, if L>L*, 



(9) 
(10) 



where L* = £oN T /T ^> £n- In a more general case of 
weak links with N conducting channels in the expres- 
sion for L* one should set 



T 



£r„. 



In what follows we will consider the most interest- 
ing limit N < N r and L<L*. Due to Eq. (9) in this 
case the dependence I e / (& x ) is fully determined by 
the current-phase relation for QPS which can be found 
by means of Eq. (6) with/ //6 = -(2e/h)dQ l fj b /d(p.We 
also note that quantum fluctuations of ip can easily be 
suppressed, e.g. by choosing the ring inductance suffi- 
ciently small. In accordance with the charge-phase un- 
certainty relation the number of Cooper pairs passing 
through the junction should strongly fluctuate thereby 
suppressing charging effects in the contact area. 

Let us now directly evaluate the supercurrent in 
both even and odd ensembles. As it was demonstrated 
above, it is sufficient to find the currents If/ b along with 
the difference of thermodynamic potentials f2j — flf = 
Qbf and combine the results with Eq. (6). The value 
If/b can be determined by means of a general formula 
[151 
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and Qbf is evaluated as a sum of the ring and the 
weak link contributions, Qb / 



fijT is found by integrating If/ b {tp) over the phase 
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ip, while the ring contribution III J a l l° w en °ugh T 



bf 

is defined by the standard expression [2,4,5] fiW ~ 
j,yA 1/2 T 3/2 exp(-A/T), where v is the density of 
states at the Fermi level, V = sL is the ring volume 
and A is the superconducting energy gap. 

3.2. Rings with quantum point contacts 

Let us first consider the limit of short QPC with d — > 

0. In this case one has Wi{u) = (2/T,)(l + uj 2 /A 2 ) - 

1, where % is the transmission of the i-th channel. 
Combining this formula with Eq. (11), after a simple 
calculation from Eq. (6) one finds [7] 
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3.3. Rings with transparent SNS junctions 
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where Ei(ip) = Ay/ 1 — % sai 2 '(ip/2). The term in the 
square brackets is specific for canonical ensembles and 
should be substituted by one for grand canonical ones. 
For N = 1 and at T = this term reduces to unity 
for even ensembles and to zero for odd ones. In other 
words, in the latter case PC will be completely blocked 
by the odd electron and no current will flow in the 
system despite the presence of a non-zero external flux 
threading the ring. In QPC with several conducting 
channels the current through the most transparent one 
will be blocked by the odd electron and, hence, the 
effect can remain significant also for QPC with N > 1. 

The physics of this blocking effect was discussed in 
details in Ref. 7. It is related to a complete cancellation 
of the contributions to the supercurrent coming from a 
pair of discrete Andreev levels with energies E±(ip) = 
±£i (<£>). Indeed, for the even ensemble at T = only 
the lower Andreev level E- (ip) is occupied, while the 
upper one E- (ip) remains empty. Accordingly, only the 
lower Andreev level contributes to the supercurrent in 
this case, and the result coincides with the standard 
grand canonical expression [16,17,18,19]. In contrast, 
for the odd ensemble both levels E±(ip) are occupied 
at T — and contribute to the supercurrent in equal 
measure. These contributions, however, enter with op- 
posite signs and completely cancel each other thus re- 
ducing the net supercurrent down to zero. We also note 
that a similar result was recently obtained in Ref. [20] 

It is important to emphasize that the above physi- 
cal picture is based on the assumption that the super- 
current is determined only by the contributions of two 
discrete Andreev levels E±(ip). This is the case pro- 
vided (a) quasiparticle states with energies above A 
do not contribute to the supercurrent and (b) there 
are no more discrete Andreev levels inside the junc- 
tion. Both these conditions are met only for symmet- 
ric and extremely short QPC with d — ► 0. In a general 
case d ^ at least one of these conditions is violated, 
and no exact compensation of the supercurrent by the 
odd electron countercurrent is anymore possible. The 
physics then becomes even more interesting, as it will 
be demonstrated below. 



From now on let us lift the condition d — > 0, i.e. 
consider SNS junctions with a non-zero thickness of 
the normal metal. Since the parity effect becomes more 
pronounced with increasing transmission of the sys- 
tem, let us first restrict our analysis to fully transpar- 
ent SNS junctions. In this case the function Wi(u) = 
W(uo) is the same for all channels. It reads 

„„ » [2u 2 \ , (2L0d\ 



2u) I u) 2 . , ( 2ujd\ 



(13) 



Substituting this function into (11) and repeating the 
the whole calculation as above, in the limit T — > we 
obtain 
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for short SNS junctions y = dA/hvF <C 1 and 
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(14) 
(15) 

nd r ' nd V N J (16) 

for long ones d>^o~ Hvf/A. These results apply for 
— n < ip < n and should be 27r-periodically continued 
otherwise. The term containing ln(l/y) in Eq. (15) for 
7 e is written with the logarithmic accuracy and is valid 
for ip not too close to ip = ±7T. 

At T — the current I e again coincides with that 
for the grand canonical ensembles [21,22]. At the same 
time for odd ensembles we observe no blocking effect 
anymore, but rather a non-zero current I for any non- 
zero value of d, cf. Eqs. (15-16). As before, at T = 
the countercurrent produced by the odd electron ex- 
actly compensates the current of the next lower An- 
dreev level. However, the contribution of continuous 
spectrum (in the case d <C Co) or of continuous spec- 
trum and all other discrete Andreev levels (in the case 
d 3> Co) remains uncompensated and is responsible for 
the non-zero current I (15-16). 

It is important to emphasize that such a current 
causes a jump on the current-phase dependence at ip = 
0, and the direction of I at sufficiently small ip is al- 
ways opposite to that of I e - This property implies that 
in the case of odd ensembles the minimum energy (zero 
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current) state occurs not at ip = 0, but at some other 
value of the phase difference. For instance, for d>^o 
the "saw tooth" current-phase relation is shifted by the 
value tt/N and, hence, the minimum Josephson energy 
(zero current) state is reached at ip = ±ir/N. One con- 
cludes that the "7r/iV-junction" state is realized in this 
case. This non-trivial state has a number of specific 
features. For instance, in the particular case N = 2 the 
current-phase relation I {<p) turns 7r-periodic. In addi- 
tion, for any N > 1 the ground state of the system ip — 
±tt/N is a twofold degenerate one within the interval 
— 7r < ip < 7r. Here we note that a similar behavior is 
expected for SNS junctions formed by d-wave super- 
conductors [24]. 

Let us also recall that the 7r-junction state can be 
realized in SNS structures by driving the electron dis- 
tribution function in the contact area out of equilib- 
rium [25,26,27]. Here, in contrast, the situation of a 7T- 
or 7r/iV-junction is achieved in thermodynamic equilib- 
rium. Despite this drastic difference, there also exists 
a certain physical similarity between the effects dis- 
cussed here and in Refs. [25,26,27]: In both cases the 
electron distribution function in the weak link deviates 
substantially from the Fermi function. It is this devi- 
ation which is responsible for the appearance of the 
7r-junction state in both physical situations. 

The current-phase relations I e (tp) and I (ip) for ar- 
bitrary values of the parameter y can be computed nu- 
merically. Examples are presented in Fig. 2. In the case 
of odd ensembles one clearly observes the current jump 
at ip — 0. As it is obvious from Eq. (15), this feature 
disappears only in the QPC limit y — » 0. 



3.4. Spontaneous currents in the ground state 

Perhaps the most interesting property of super- 
conducting rings interrupted by a 7r-junction is the 
possibility to develop spontaneous supercurrent in the 
ground state [23]. Although this feature is inherent to 
any type of 7r-junctions, in the case of the standard 
sinusoidal current-phase relation such spontaneous 
supercurrents can occur only for sufficiently large val- 
ues of the ring inductance C [23]. In contrast, in the 
situation studied here the spontaneous current state is 
realized for any inductance of the ring because of the 
non-sinusoidal dependence I {<p) (16) [28]. 

In order to demonstrate that let us assume that no 
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Fig. 2. The zero temperature current-phase relations I e {<f) 
and I (<p) {— t < ip < I") for N = 1 and different values of 
the parameter y = dA / hvp . 

external flux is applied to our system. Then at T — > 
the energy of an SNS ring with odd number of elec- 
trons can be written in the form 



(17) 



where $ is the flux related to the circular current flow- 
ing in the ring. Minimizing this energy with respect 
to $, one easily observes that a non-zero spontaneous 
current 



(18) 



should flow in the ground state of our system. This 
is yet one more remarkable consequence of the parity 
effect: Just by changing M from even to odd one can 
induce non-zero PC without any external flux & x . In 
the limit of small inductances C <S ^od/evpN - which 
is easy to reach in the systems under consideration - the 
value of / does not depend on the number of channels 
N and is given by the universal expressions can easily 
be derived from Eqs. (15-16): 



(19) 
(20) 



For intermediate values of the parameter y the ampli- 
tude of the current I sp can be evaluated numerically. 
The results are displayed in Fig. 3. In agreement with 
Eq. (19) the spontaneous current I sp increases linearly 
with d at small d, reaches its maximum value I ma x ~ 
0AeA/h at d ~ £o and then decreases with further 
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Fig. 3. The spontaneous current amplitude I sp as a function 
of the parameter y = dA/hvp at T = 0. In the inset, the 
same function is shown on the log — log scale. Dashed lines 
indicate the asymptotic behavior of I S p{y) in the limits of 
small and large y. 

increase of d approaching the dependence (20) in the 
limit of large d. For typical BCS superconductors this 
maximum current can be estimated as J max ~ 10 4- 100 
nA. Spontaneous currents of such a large magnitude 
can be directly detected in modern experiments. 



3.5. Junctions with resonant transmission 

As we have already emphasized above, the effects dis- 
cussed here are mostly pronounced for SNS junctions 
with few conducting channels N ~ 1. The most promis- 
ing candidates for practical realization of such struc- 
tures are junctions with iV-layers formed by carbon 
nanotubes or organic molecules. In both these cases it 
is not always easy to achieve good quality contacts be- 
tween nanotubes/molecules and superconducting elec- 
trodes. In other words, although the electron transport 
in the JV-metal itself can be ballistic, electron scatter- 
ing at iVS'-interfaces may be substantial and may sig- 
nificantly modify the above results. 

At the first glance, in the limit of small transmis- 
sion of iVS-interfaces the influence of the parity effect 
on PC should become weak and, hence, can be disre- 
garded. Below we will demonstrate that this conclu- 
sion - although correct under certain conditions - does 
not necessarily apply to nanojunctions for which size 
quantization of energy levels is an important effect. 
This effect gives rize to resonant transmission of elec- 



trons across the junction and considerably modifies the 
physical situation. 

To proceed we will adopt the standard model of an 
SINIS junction ("7" stands for the insulating layer) 
and assume that transmissions of both NS interfaces 
are much smaller than unity. In this limit electron 
transport across the junction is mainly due to resonant 
tunneling through discrete energy levels inside the nor- 
mal metal. For simplicity we will restrict our analysis 
to a single channel junction N = 1. If needed, gen- 
eralization of our analysis to the case N > 1 can be 
performed in a straightforward manner. 

As before, we will make use of the general expression 
(11) where the function W(ui) now takes the form [15] 



W{uj) = — — — — cos(2fc F d+, 



(21) 
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Here 7i,2 and 72-i,2 = 1 — 71,2 are respectively trans- 
mission and reflection coefficients of the two NS inter- 
faces, f2 = VA 2 + lo 2 and 2kpd + <f> is the phase accu- 
mulated by electrons during the motion between the 
barriers. 

In accordance with our assumption below we shall 
consider the limit 7i,2 <C 1. The most interesting phys- 
ical situation is realized in the limit of short SINIS 
junction d <g £0 ~ Ruf/A. Under the above conditions 
the size energy level quantization in the normal metal 
becomes an important effect. In the case of a one di- 
mensional metal of length d, the level spacing in the 
vicinity of the Fermi energy is 8e ~ hvp/d. Hence, the 
condition for the short junction regime can also be rep- 
resented in the form Se 3> A. Tunneling of electrons 
into the reservoirs causes a non-zero linewidth of the 
energy levels which is proportional to Ti^Se. This value 
is much smaller than 8e, hence, the resonances remain 
well separated. In this situation it suffices to take into 
account only the closest to the Fermi energy level in- 
side the normal metal. 

Let us introduce the energy of the resonant level tn 
by means of the relation cos(2fcF+0) ~ —1 + 1/2 (|&) , 
where Se — hvp/2d, and the tunneling rates Ti^/fo = 
7i,2<5e//t. Expanding the cosh and sinh terms in Eq. 
(21) to the leading order in A/Se and substituting the 
result in Eq. (11), we obtain 
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a b 




Fig. 4. Andreev levels in a single mode SNS junction with 
d = hv F /A: (a) Ti, 2 = 1 and (b) T M « 1 and A/T = 0.5. 
In both cases T = 1. 

r = £ T V- A 2 Tsiny /ooN 

//6 ft ^ e 2 (v5)+^ /6 (l + 40T/T mM )' 1 j 
where 

e(p) = Ax/l-Tsin 2 (^/2), T = Ti + T 2 , T max = 
4riT2/r 2 and the total transmission probability at the 
Fermi energy T is given by the Breit-Wigner formula 

T _ rxr 2 
(«) 2 + |r 2 

It is interesting to point out that the expressions de- 
fined by Eqs. (22,23) can also be derived from the An- 
derson model for resonant tunneling through a single 
impurity center between two superconductors in the 
limit of zero on-site interaction. 

It follows from Eqs. (22, 23) that - although the 
transparencies of both barriers are low - the total 
transmission T and, hence, the Josephson current 
show sharp peaks provided the Fermi energy becomes 
close to the bound states inside the junction. On the 
other hand, Eqs. (22, 23) demonstrate that even in the 
vicinity of resonances the behavior of the Josephson 
current as a function of the phase difference <p and 
temperature T can substantially deviate from that for 
transparent SNS junctions. 

In order to understand the physical reasons for such 
a difference it is instructive to compare the structure 



of discrete Andreev levels for ballistic (71,2 = 1) SNS 
junctions with that for SINIS junctions with 71 j2 <C 
1, see Fig. 4. The spectrum of an SINIS system con- 
sists only of a single non-degenerate state £o(v>) in the 
interval < £o < A (Fig. 4b). As a result, the behav- 
ior of eo(v) at small ip is smooth and the derivative of 
so with respect to ip has no jump at if = 0. In contrast, 
in the case of ballistic SNS junctions discrete levels 
become split at arbitrary small values of the phase ip 
(Fig. 4a) . Hence, the derivative of the lowest Andreev 
level with respect to ip acquires a jump at ip = 0. As 
we have already discussed, this feature is crucial for 
the presence of spontaneous currents in superconduct- 
ing rings with arbitrary small inductances C and with 
odd number of electrons. Since this feature is absent in 
the case of SINIS junctions the spontaneous current 
in the ground state of the system can only develop at 
not very small ring inductances. 

Analytical expressions for the energy of the subgap 
state (obtained from the poles of the summand in Eq. 
(22)) and for the parity-dependent currents I e / can 
be derived in various asymptotic regimes. At a given 
value of the total transmission T one can distinguish 
two limiting cases: A/r « 1 and A/T ^> 1. In the 
case of a wide resonance A/F <C 1, the energy of the 
subgap state is Eo(ip) = Ayl — Tsin 2 (ip/2) and I e / 
are given by the expressions derived for a single mode 
QPC with an effective transmission T. In the opposite 
narrow resonance limit A/F> 1 we get 

so iy) = ^«) 2 + ilV 1 ^ T sin2 fa/ 2 ) • 

The currents I e / a are again given by the expressions 
obtained for QPC with the substitution of A by 
\J (ffl) 2 + jT 2 . Note that in both these asymptotic 
regimes only the discrete spectrum contributes to the 
currents I e /o- Thus, for these regimes the Josephson 
current will be completely blocked by the odd electron 
in the limit of low temperatures. 

At intermediate values of the parameter A/T the 
continuous spectrum will contribute to the Josephson 
current. The current-phase relations I e / (ip) can be 
easily evaluated numerically for arbitrary parameter 
values. The results of these calculations - partially rep- 
resented in Fig. 5 - clearly demonstrate that at suffi- 
ciently low temperatures the "7r-junction" state should 
be realized in the case of odd number of electrons. Ac- 
cording to our expectations, however, in the case of 
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Fig. 5. Zero temperature current-phase relations I e / {<p) 
for T = 0.9, Ti = T2 = To and different values of the 
parameter A/Fo. 

resonant SINIS junctions the current has no jump at 
if = and the behavior of I {f) is more similar to that 
of "conventional" ^-junctions [23]. 

Let us emphasize again that all the above results for 
SINIS junctions with small interface transmissions 
were derived in the limit of short junctions d <g £o- Res- 
onant effects are also important in the opposite limit 
d >• £o [15]. However, in the latter case even exactly at 
resonance the Josephson current remains small being 
proportional to the smallest of the two transmissions 
Ti and T2. Hence, the parity effect in the limit d 3> £0 
will also be small in equal measure. For this reason the 
case of long junctions with resonant transmission will 
not be considered here. 



4. Conclusions 

In this paper we have demonstrated that at suffi- 
ciently low temperatures superconducting parity ef- 
fect may strongly influence equilibrium persistent cur- 
rents in isolated superconducting nanorings contain- 
ing a weak link with few conducting modes. An odd 
electron, being added to the ring, occupies the lowest 
available Andreev state and produces a countercurrent 
circulating inside the ring. For a single channel quan- 
tum point contact at T = this countercurrent exactly 
compensates the supercurrent I e produced by all other 
electrons and, hence, yields complete blocking of PC for 
any value of the external magnetic flux. In supercon- 



ducting nanorings with embedded normal metal the 
odd electron countercurrent can "overcompensate" I e 
and a novel "it /N-j unction" state occurs in the system. 
Changing the electron parity number from even to odd 
results in spontaneous supercurrent in the ground state 
of such rings without any externally applied magnetic 
flux. This and other novel effects predicted here can 
be directly tested in modern experiments. They can 
also be used for practical realization of superconduct- 
ing qubits. 



Acknowledgments This work is part of the Kom- 
petenznetz "Funktionelle Nanostructuren" supported 
by the Landestiftung Baden- Wurttemberg gGmbH. 
This work has also been supported by the European 
Community's Framework Programme NMP4-CT- 
2003-505457 ULTRA- ID "Experimental and theoreti- 
cal investigation of electron transport in ultra-narrow 
1-dimensional nanostructures" . 



References 

[I] D.V. Averin and Yu.V. Nazarov, Phys. Rev. Lett. 69 
(1992) 1993. 

[2] M.T. Tuominen, J.M. Hergenrother, T.S. Tighc, and 
M. Tinkham, Phys. Rev. Lett. 69 (1992) 1997. 

[3] P. Lafargc, P. Joyez, D. Estcvc, C. Urbina, and M.H. 
Devoret, Phys. Rev. Lett. 70 (1993) 994. 

[4] B. Janko, A. Smith, and V. Ambegaokar, Phys. Rev. 
B 50 (1994) 1152. 

[5] D.S. Golubcv and A.D. Zaikin, Phys. Lett. A195 (1994) 
380. 

[6] D.V. Averin and Yu.V. Nazarov, Physica B203 (1994) 
310. 

[7] S.V. Sharov and A.D. Zaikin, Phys. Rev. B 71 (2005) 
014518. 

[8] K. Rang, Europhys. Lett. 51 (2000) 181. 

[9] H.-J. Kwon and V.M. Yakovcnko, Phys. Rev. Lett. 89 
(2002) 017002. 

[10] A.D. Zaikin, D.S. Golubcv, A. van Ottcrlo, and G.T. 
Zimanyi, Phys. Rev. Lett. 78 (1997) 1552. 

[II] D.S. Golubcv and A.D. Zaikin, Phys. Rev. B 64 (2001) 
014504. 



8 



[12] A. Bezryadin, C.N. Lau, and M. Tinkham, Nature 
(London) 404 (2000) 971. 

[13] C.N. Lau, N. Markovic, M. Bockrath, A. Bezryadin, 
and M. Tinkham, Phys. Rev. Lett. 87 (2001) 217003. 

[14] The relation (7) is strictly applicable only for grand 
canonical ensembles. However, at this point the 
difference between canonical and grand canonical 
ensembles is unimportant and can be disregarded. 

[15] A.V. Galaktionov and A.D. Zaikin, Phys. Rev. B 65 
(2002) 184507. 

[16] I.O. Kulik and A.N. Omel'yanchuk, Sov. J. Low Temp. 
Phys. 4 (1978) 142. 

[17] W. Haberkorn, H. Knauer, and J. Richter, Phys. Stat. 
Solidi (A) 47 (1978) K161. 

[18] A. Furusaki and M. Tsukada, Physica B165-166 (1990) 
967. 

[19] C.W.J. Beenakker and H. van Houten, Phys. Rev. Lett. 
66 (1991) 3056. 

[20] M. Hayashi and H. Ebisawa, Phys. Rev. B 67 (2003) 
014524. 

[21] I.O. Kulik, Sov. Phys. JETP 30 (1970) 944. 

[22] C. Ishii, Progr. Thcor. Phys. 44 (1970) 1525. 

[23] L.N. Bulacvskii, V.V. Kuzii, and A. A. Sobyanin, JETP 
Lett. 25 (1977) 290. 

[24] Yu.S. Barash, A.V. Galaktionov and A.D. Zaikin, 
Phys. Rev. B 52 (1995) 665. 

[25] A.F. Volkov, Phys. Rev. Lett. 74 (1995) 4730. 

[26] F.K. Wilhelm, G. Schon, and A.D. Zaikin, Phys. Rev. 
Lett. 81 (1998) 1682. 

[27] S.K. Yip, Phys. Rev. 58 (1998) 5803. 

[28] For small C — > non-vanishing spontaneous 
currents in the ground state may as well exist in 
systems containing SNS junctions formed by d-wave 
superconductors [24]. Also in that case this feature is 
due to the non-sinusoidal current-phase relation in such 
junctions. 



